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1. Introduction 

In recent years, the mixed-spin Ising models belong to the most actively studied lattice-statistical models in 
the statistical and solid-state physics as they often exhibit unpredictable and rather complex critical behavior. 
In particular, the Ising systems containing spins of different magnitudes are still relatively simple but useful 
models that bring a deeper insight into the magnetic behavior of certain ferrimagnetic materials, which are of 
great technological importance due to their possible applications in thermomagnetic recording [1]. With regard 
to this, the investigation of ferrimagnetic behavior of the mixed-spin Ising models has become a very active 
research field over the last few decades. Despite the intensive studies, however, there are only few examples of 
exactly solved mixed-spin Ising models, yet. Using a generalized form of the decoration-iteration and star-triangie 
transformations, the mixed spin-1/2 and spin-5 (5* > 1) Ising models on the honeycomb, diced and decorated 
honeycomb lattices were exactly been treated by Fisher [2] and Yamada [3] many years ago. These mapping 
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transformations were later on further generalized in order to account also for the single-ion anisotropy terms. 
Actually, the procedure based on generalized mapping transformations was recently employed to obtain exact 
results for the mixed-spin Ising models with the uniaxial single-ion anisotropy on the honeycomb lattice [4-6], 
bathroom-tile (4-8) lattice [7], as well as, several decorated planar lattices [8-10]. On the other hand, the mixed- 
spin Ising models that account also for the biaxial single-ion anisotropy have exactly been solved just on the 
honeycomb [11] and diced [12] lattices, so far. It is worthy to mention that main difficulties, which emerge when 
treating the Ising models refined by the biaxial single-ion anisotropy term, originate from a presence of the x- 
and j/-components of spin operators introducing to a spin system quantum fluctuations. 

Owing to these facts, the Ising models accounting for the biaxial single-ion anisotropy have recently been a subject 
matter of many approximative theories as well. The ground state of the general spin-S Ising model with the 
biaxial single-ion anisotropy was investigated by establishing an effective mapping to the transverse Ising model 
[13] and by using the effective-field theory with self-spin correlations [14]. It is worthwhile to remark that the 
finite-temperature properties of the single-spin Ising models with the biaxial single-ion anisotropy have also been 
investigated within the framework of the variational mean-field treatment [15], the random phase approximation 
[16, 17], the linked-cluster series expansion [18] and the standard effective-field theory with correlations based 
on the differential operator [19] or probability distribution [20] techniques. Contrary to this, the effect of biaxial 
single-ion anisotropy on magnetic properties of the mixed-spin Ising models was less frequently studied and thus, 
it is still not fully understood, yet. Except few exactly solved cases mentioned earlier [11, 12], the mixed-spin Ising 
models with the biaxial single-ion anisotropy term have been explored just within the conventional effective-field 
theory based on the differential operator [21] and probability distribution [22] techniques. 

With all this in mind, in this paper we will investigate the magnetic properties of the mixed-spin Ising model 
on the Bethe lattice with both uniaxial and biaxial single-ion anisotropy terms. Exact solution for this model 
system will be obtained by combining two accurate mapping transformations with the exact method based on 
recursion relations. First, the star-triangle transformation is used to connect the mixed-spin Ising model on the 
Bethe lattice with the coordination number g=3 to its equivalent spin-1/2 Ising model on the triangular Husimi 
lattice. Next, the spin-1/2 Ising model on the triangular Husimi lattice is subsequently mapped by the use of 
triangle-star transformation to the simple spin-1/2 Ising model on the Bethe lattice. It is well known that this 
latter model can be exactly treated using exact recursion relations [23], which will help us to complete our exact 
calculation for the original mixed-spin Ising model on the Bethe lattice. 

The outline of this paper is as follows. In Section 2, basic steps of the exact mapping transformations will be 
explained together with some details, which enable to find several exact results for the mixed-spin Ising model 
on the Bethe lattice with both uniaxial as well as biaxial single-ion anisotropy terms. In Section 3, the most 
interesting results will be presented and discussed for ground-state and finite-temperature phase diagrams, the 
total and sublattice magnetization. In this section, our attention will also be focused on a possibility of observing 
compensation phenomenon. Finally, some concluding remarks are drawn in Section 4. 
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Figure 1. The fi gurc on the left shows the mixed-spin Bethe lattice with the coordination number g— 3, which consists of two 
inequivalent interpenetrating sublattices A and B constituted by the spin-1/2 (empty circles) and spin-1 (filled circles) 
atoms, respectively. The mixed-spin Ising model on the Bethc lattice is mapped via the star-triangle transformation on 
the spin-1/2 Ising model on the triangular Husimi lattice depicted in the central part of this figure. Subsequently, the 
spin-1/2 Ising model on the triangular Husimi lattice is mapped through the triangle-star transformation on the simple 
spin-1/2 Ising model on the Bethe lattice, which is formally divided into two equivalent interpenetrating sublattices A 
and B whose sites are shown as empty and gray circles, respectively. 

2. Model and its exact solution 

Let us begin by considering a mixed spin-1/2 and spin-1 Ising model on the Bethe lattice with the coordination 
number q—3, which is schematically illustrated on the left-hand-side of Fig. 1. As shown in this figure, the 
mixed-spin Bethe lattice contains a central spin-1 site which has q nearest-neighbor spin-1/2 sites forming the 
first generation. The second generation is formed by connecting g-1 new spin-1 sites to each site of the first 
generation and so on. Thus, the mixed-spin Bethe lattice consists of two inequivalent interpenetrating sublattices 
A and B, which are occupied by the spin-1/2 and spin-1 atoms depicted as empty and filled circles, respectively. 
Assuming the Ising-type exchange interaction J between nearest-neighboring spins from two different sublattices 
A and B, the total Hamiltonian of the investigated spin system is given by 

2iV AT AT 

(i.j) i=i i=i 

where A/ is a total number of lattice sites at each sublattice, (a = x, y, z) and it| denote the spatial components 
of the usual spin-1 and spin-1/2 operators, respectively. The first summation is carried out over all nearest- 
neighboring spin pairs, while the other two summations run over all sites of the sublattice B. The last two 
terms D and E, which affect just the spin-1 atoms from the sublattice B, measure a strength of the uniaxial and 
biaxial single-ion anisotropy, respectively. It should be mentioned here that by neglecting the biaxial single-ion 
anisotropy, i.e. setting iJ=0 in Eq. (1), our model reduces to the mixed-spin Ising model with the uniaxial 
single-ion anisotropy exactly solved in several previous works [24-26]. In this work, we will therefore mainly 
investigate the effect of the biaxial single-ion anisotropy, which infiuences magnetic properties of the model under 
consideration in a crucial manner. Really, the E term related to the biaxial single-ion anisotropy might cause 
non-trivial quantum effects since it introduces the x and y components of spin operators into the Hamiltonian 
(1) and thus, it is responsible for the onset of local quantum fiuctuations that are missing in the Ising models 
accounting for the uniaxial single-ion anisotropy only. 
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Now, let US turn our attention to main points of the mapping transformation method, which enables an exact 
treatment of the model under investigation. By introducing site Hamiltonians, the total Hamiltonian (1) can be 
written as a sum over all site Hamiltonians 

N 

'H = J2'Hi, (2) 

i = l 

where each site Hamiltonian Hi involves all the interaction terms associated with the spin-1 atom residing on the 
ith site of the sublattice B 

m = -StE, - {StfD - [{S^f - (5f )^]£;, (3) 

with Ei = J((Tfi + (7^2 + o'a)- Because the Hamiltonians (3) at different sites commute, i.e. \Hi,'Hj\ = is valid 
for each i ^ j, the partition function of this spin system can be partially factorized and consequently rewritten 
in the form 

N 

^ = EnTrs.exp(-/?-«0, (4) 

{a} i = l 

where (3 = l/(fcsr), ks is Boltzmann's constant, T is the absolute temperature, X]{o-} denotes a summation 
over all possible spin configurations on the sublattice A and Trg^ stands for a trace over spin degrees of freedom 
of ith spin from the sublattice B. It should be mentioned that a crucial step of our exact procedure represents 
calculation of the expression Trs^ exp(— /37ii). For this purpose, it is useful to write the site Hamiltonian (3) in 
the matrix form 



v 



-Ei- D Q E 

E E^-D 



(5) 



in a standard basis of functions |±1), |0) corresponding, respectively, to the three possible spin states Sf — ±1, of 
ith atom from the sublattice B. After a straightforward diagonalization of the site Hamiltonian (5), the expression 
for a partial trace over the spin states of ith spin immediately implies a possibility of performing the generalized 
star-triangle mapping transformation 



Trs. exp(-/3-W0 = 1 -f 2exp(/3D) cosh(^/3y^J2(afi + af^ + + E^^ 
= A exp |^/3i?((jfiO-f2 + f^i2cri3 + crfgcra) j , 



(6) 



which replaces the partition function of a star, i.e. the four-spin cluster consisting of one central spin-1 site and 
its three enclosing spin-1/2 sites, by the partition function of a triangle, i.e. the three-spin cluster comprising of 
three spin-1/2 sites to be located in the corners of equilateral triangle (see Fig. 1). The physical meaning of the 
mapping (6) is to remove all interaction parameters associated with the central spin-1 atom and to replace them 
by an effective interaction R between the outer spin-1/2 atoms. It is noteworthy that both mapping parameters 
A and R are "self-consistently" given directly by the transformation (6), which must hold for any combination 
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of the spin states of three spin- 1/2 atoms. After taking into account all possible spin configurations of the three 
spin-1/2 atoms, one obtains from the star-triangle transformation (6) just two independent equations that directly 
determine so far not specified mapping parameters A and R 

A=(v,Viy^\ (3R = ln(^). (7) 

In above, we have defined the functions Vi and V2 in order to express the transformation parameters A and R in 
a more simple and elegant form: 

Vi = 1 + 2 exp(/?D) cosh (^/3y^(3J/2)2 +i;2^ , 

V2 = l + 2exp{PD) cosh(^/?y^(J/2)2 -f £2^. (8) 

When the star-triangle mapping (6) is performed at each site of the sublattice B, the original mixed-spin Ising 
model on the Bethe lattice is mapped on the spin-1/2 Ising model on triangular Husimi lattice with the effective 
nearest-neighbor interaction R given by the "self-consistency" condition (7)-(8) (see Fig. 1). As a matter of fact, 
the substitution of the mapping transformation (6) into the partition function (4) establishes the relationship 

Z{(3, J, D, E) = A'^2Husimi(/3, R), (9) 

between the partition function Z of the mixed-spin Ising model on the three-coordinated Bethe lattice and 
respectively, the partition function ^Husimi of the corresponding spin-1/2 Ising model on the triangular Husimi 
lattice. By the triangular Husimi lattice is meant a deep interior of the Husimi tree [27], which is built up 
from corner-sharing triangles attached to each other by their vertices as it is schematically shown in the central 
part of Fig. 1. Even though the spin-1/2 Ising model on the triangular Husimi lattice can already be exactly 
solved with the help of exact recursion relations [28-30] , it seems for us more useful to perform another mapping 
transformation that relates the investigated model system to the simple spin-1/2 Ising model on the Bethe lattice 
for which few exact analytical results are available in the literature [23]. 

In the next step, let us therefore perform the triangle-star transformation by inserting the spin-1/2 atom into 
each triangle of the triangular Husimi lattice (see Fig. 1). This procedure will allow us to express the partition 
function of the spin-1/2 Ising model on the triangular Husimi lattice in terms of the partition function of the 
simple spin-1/2 Ising model on the Bethe lattice. However, it is more advisable to start from the Hamiltonian 
of the simple spin-1/2 Ising model on the three-coordinated Bethe lattice and to show that its partition function 
can really be connected to the partition function of the spin-1/2 Ising model on the triangular Husimi lattice 
after performing the appropriate star-triangle transformation. The spin-1/2 Ising model on the Bethe lattice 
schematically shown on the right-hand-side of Fig. 1 can be defined through the Hamiltonian 



(10) 
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where we have formally distinguished the Ising spins erf = ±1/2 and jij = ±1/2 to be shown as empty and gray 
circles, respectively, in order to divide the three-coordinated Bethe lattice into two equivalent interpenetrating 
sublattices A and B. Then, the procedure that leads to the star-triangle mapping transformation can be repeated 
once again. The total Hamiltonian (10) can be firstly rewritten as a sum of A'' site Hamiltonians 



where each site Hamiltonian Hi 



(11) 



(12) 



involves all the interaction terms of jth spin from the sublattice B. With regard to the above definition of site 
Hamiltonians (12), it is possible to write the following expression for the partition function ^Bothc of the spin-1/2 
Ising model on the Bethe lattice 

JV 

~ " (13) 



{rr} i = l fjj=±l/2 



where the symbol X]{o-} denotes a summation over all possible spin configurations on the sublattice A and the 
second summation is carried out over particular spin states of jth spin from the sublattice B. When this 

particular summation is explicitly carried out, one gains the expression that in turn implies a possibility of 
performing the usual star-triangle mapping transformation [2] 



^ exp{-PHj) = 2 cosh (o-|i ± a^a ± (Tjs)^ = B exp [PR (cr|io-|2 + '^h'^h + '^h'^ji)] 



(14) 



Mj=±l/2 



with so far not specified mapping parameters B and R. The self-consistency condition yields for the mapping 
parameters B and R the following relations 



B = 2 tl cosh ( ^l3Kj cosh'* (jPK ) , 



/3R = In 



2cosh(i/3A') 



(15) 



When the mapping transformation (14) with appropriately chosen mapping parameters (15) is performed at each 
site of the sublattice B, the spin-1/2 Ising model on the three-coordinated Bethe lattice with the nearest-neighbor 
interaction K is mapped onto the spin-1/2 Ising model on the triangular Husimi lattice with the effective nearest- 
neighbor interaction R. By substituting the transformation formula (14) into the partition function (13), one 
indeed obtains a rather simple mapping relationship between the partition functions of the spin-1/2 Ising model 
on the three-coordinated Bethe lattice and respectively, the spin-1/2 Ising model on triangular Husimi lattice 



ZBctha{P,K) = B^ ZuusimiiP , R) ■ 



(16) 



Of course, it is also possible to find inverse triangle-star transformation that expresses the partition function 
of the spin-1/2 Ising model on triangular Husimi lattice in terms of that one for the corresponding spin-1/2 
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Ising model on the three-coordinated Bethe lattice, namely, Eq. (16) can simply be inverted. Consequently, the 
partition function of the mixed spin- 1/2 and spin-1 Ising model on the Bethe lattice can also be straightforwardly 
calculated from the partition function of the corresponding spin-1/2 Ising model on the Bethe lattice. After 
eliminating the partition function of the spin-1/2 Ising model on triangular Husimi lattice from Eqs. (9) and (16), 
one actually obtains a simple mapping relationship 



which relates the partition function of the mixed spin-1/2 and spin-1 Ising model on the Bethe lattice to that one 
of the corresponding spin-1/2 Ising model on the Bethe lattice. Besides, the effective nearest-neighbor interaction 
of the corresponding spin-1/2 Ising model on the three-coordinated Bethe lattice can readily be obtained by 
eliminating the effective interaction (3R from Eqs. (7) and (15). This procedure enables to express the mapping 
parameter K, which determines a strength of the nearest-neighbor interaction of the corresponding spin-1/2 Ising 
model on the Bethe lattice, solely as a function of the temperature and interaction parameters J, D and E 
originally included in the Hamiltonian (1) 



The sign ambiguity to emerge in the previous equation reflects the zero-field invariance of the spin-1/2 Ising 
model on the Bethe lattice with respect to the transformation K — > ~K, because the relevant sign change does 
not affect a strength (absolute value) of the effective coupling f3K. As a result, the mixed spin-1/2 and spin- 
1 Ising model on the Bethe lattice can alternatively be mapped either to the ferromagnetic (plus sign) or the 
antiferromagnetic (minus sign) spin-1/2 Ising model on the Bethe lattice, which both have the identical critical 
temperature /3c = 2 In 3 for the same strength of the effective interation (from here onward, we will only utilize 
the solution with plus sign for simplicity). It is worthwhile to remark, moreover, that the above equation in fact 
completes our exact calculation, since the partition function of the mixed spin-1/2 and spin-1 Ising model on 
the Bethe lattice with both uniaxial as well as biaxial single-ion anisotropy terms can be now extracted from the 
well-known exact result for the partition function of the corresponding spin-1/2 Ising model on the Bethe lattice 
[23] unambiguously given by the effective nearest-neighbor interaction (18). 

Exact results for phase diagrams and other relevant physical quantities now follow straightforwardly. For instance, 
it can be easily understood from the mapping relation (17) that the mixed-spin- Ising model on the Bethe lattice 
becomes critical if and only if the corresponding spin-1/2 Ising model on the Bethe lattice becomes critical as well, 
since the mapping parameters A and B given by Eqs. (7)-(8) and (15) are analytic function in the whole region 
of the interaction parameters. From this point of view, it is sufficient to compare the effective nearest-neighbor 
interaction (18) of the spin-1/2 Ising model on the Bethe lattice with its critical value {l3cK = 21n3) 




(17) 




(18) 




(19) 
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in order to determine a critical behavior of the model under investigation (the superscript c in the above equation 
means that the inverse critical temperature /3c = i/iksTc) enters into the expressions Vf and V2 given by 
Eq. (8) instead of P). Similarly, the mapping relations can also be used to determine the total and sublattice 
magnetizations. It can be easily proved by using the mapping theorems developed by Barry et al. [31] that the 
sublattice magnetization niA of the mixed-spin Ising model on the Bethe lattice is directly equal to the single-site 
magnetization of the corresponding spin-1/2 Ising model on the triangular Husimi lattice given by the effective 
interaction R, as well as, the single-site magnetization of the corresponding spin-1/2 Ising model on the Bethe 
lattice given by the effective interaction K 

niA = {d-t) = )Husimi,H — (^-j^ ) Bethe, K = mBetho- (20) 

In above, the symbols {...), {. . .)Husimi,fl and {. . .)Bethe,jf denote canonical ensemble average performed within 
the mixed spin-1/2 and spin-1 Ising model on the Bethe lattice, the spin-1/2 Ising model on the Husimi lattice 
with the effective interaction R and the spin-1/2 Ising model on the Bethe lattice with the effective interaction K, 
respectively. From this point of view, the sublattice magnetization uia can easily be obtained from the iteration 
procedure based on the exact recursion relations yielding 

where x is given by the stable fixed point of the recurrence relations [23] 

exp(-i^)-fexp(g|^)4_i 



exp(if )-fexp(-^):r2_, 



(22) 



On the other hand, the sublattice magnetization ms can be obtained after straightforward but a little bit cum- 
bersome calculation from the exact Callen-Suzuki identity [32] 

ms^iSt) = h''^'^'^^^-^^A, (23) 
^ \ Trs, exp(-/3H0 /' 

which enables to express the sublattice magnetization ms in terms of the formerly derived sublattice magnetization 
ruA and the triplet correlation function Ia = {d'iid-i2^i-i) 

TUB = 1^a[F{^) + F(i)] + 2tA[F{^) - 3F{^)]. (24) 

For completeness, the function F{x) is defined as follows 



^/{Jxy + £2 I + 2exp{i3D) cosh(/3-^(Ja;)2 4- E^) 
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and the triplet correlation function, which is required for a computation of the sublattice magnetization ms, 
can again be calculated within the framework of exact mapping theorems [31]. The exact mapping theorems 
give for the unknown triplet correlation function on the Bethe lattice the following result Ia = {S'iiO'i2&i3) = 
{(7fi(Tf20"i3)Husimi,fl = (o'fi '5"f2'5"i3 ) Bctho.K = tectho and they also connect it to the single-site magnetization through 
the relation 

2-3[G(|)+G(l)] 

teethe = "^Bethe ^[^^ 3 ) _ 1 )j • (26) 

In above, the newly defined function is G{x) — i tanh(^a;). 

3. Results and discussion 

Although all derivations presented in the preceding section hold for the ferromagnetic (J > 0) as well as ferri- 
magnetic (J < 0) version of the model under investigation, in what follows we shall restrict ourselves only to 
an analysis of the ferrimagnetic model with J < 0. It is worthy to mention, nevertheless, that phase diagrams 
displayed below for the ferrimagnetic model are valid without any changes also for the ferromagnetic model as a 
result of an invariance of the mapping transformations with respect to J ^ — J interchange. 
First, let us take a closer look at the ground-state behavior of the mixed spin-1/2 and spin-1 Ising model on the 
Bethe lattice. Taking into account the zero temperature limit (T ~> 0'^), one finds the following condition for a 
first-order phase transition line separating the magnetically ordered phase (OP) and the disordered phase (DP): 



The detailed analysis shows that the spin ordering appearing within OP and DP can unambiguously be defined 
through the following eigenfunctions 



|0P) = n 1 - ^) n 1 1) - sin <^ I - 1)] , (28) 

jeA ieB 

iDP)=ni±^)nio). (29) 



where the former products are carried out over all the spin-1/2 sites from the sublattice A, the latter ones run 
over all the spin-1 sites from the sublattice B and the angle determines probability amplitudes for the | + 1) and 
I — 1) spin states within OP through the relation 4> = ^ arctan (^jp^^ ■ It can easily be understood from Eq. (29) 
that all the spin-1 atoms reside within DP their "non-magnetic" spin state |0) and as a result of this, all the 
spin-1/2 atoms become completely disordered, i.e. they are completely randomly in one of their two possible spin 
states I ± 1). However, the more striking situation emerges within OP, where the spin-1 atoms are in a quantum 
superposition of two spin states j + 1) and j — 1) and simultaneously, all the spin-1/2 atoms reside their "down" 
spin state | — ^). Obviously, the biaxial single- ion anisotropy enhances the probability amplitude of the spin state 
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-1.0 




Figure 2. a) Ground-state phase diagram in the E — D plane; b) Single-site magnetizations versus the biaxial single-ion anisotropy 
E/\J\ and D/\.J\ = -1.0 at T = 0. 

I — 1) and thus, the E term effectively lowers the sublattice magnetization uib due to the quantum reduction of 
the magnetization (the sublattice magnetization m,4 is not directly affected by this term). In agreement with the 
aforedescribed ground-state analysis, the following analytical expressions for the single-site magnetization {niA, 
ma) and the total single-site magnetization m — {niA + mB)/2 follow from Eqs. (20)-(26) at the ground state: 




DP : mA = 0.0, tub = 0.0, m = 0.0. (31) 

Figs. 2 (a) and 2(b) show the ground-state phase diagram in the E-D plane and zero-temperature variations 
of the magnetization with the biaxial single-ion anisotropy when D/\J\— -1.0. It should be mentioned that by 
neglecting the biaxial single- ion anisotropy, i.e. putting E—0 into the condition (27), one recovers the phase 
boundary between OP and DP for the uniaxial single-ion anisotropy D/\J\=-1.5 that is consistent with previous 
calculations for this mixed-spin system on the Bethe lattice [5, 24-26]. In this particular case, both sublattice 
magnetizations have antiparallel orientation with respect to each other and they are also both fully saturated. It 
should be also remarked that OP corresponds in this special case to the simple ferrimagnetic phase. However, the 
ground-state behavior becomes much more complex in the presence of the non-zero biaxial single-ion anisotropy 
E. Although the sublattice magnetization niA remains at its saturation value in the whole OP, the sublattice 
magnetization ms monotonically decreases by increasing a strength of the biaxial single-ion anisotropy. Obviously, 
the decrease of the sublattice magnetization tub implies a violation of a perfect ferrimagnetic spin ordering in the 
OP due to the aforedescribed quantum reduction of the magnetization caused by the E term. 
Now, let us proceed to the finite-temperature properties of the system under investigation by considering the effect 
of uniaxial and biaxial single-ion anisotropies on the critical behavior. The dependence of the critical temperature 
Tc on the uniaxial and biaxial single-ion anisotropies is shown in Fig. 3 and Fig. 4, respectively. In both figures, 
the depicted phase transition lines separate OP from DP, which always appear for temperatures above the phase 
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transition lines. Note furthermore that the phase transitions between these two phases are of second-order with 
the standard mean-field like critical exponents. More specifically, Fig. 3 shows the critical temperature as a 
function of the uniaxial single-ion anisotropy for different values of the biaxial single-ion anisotropy. It can be 
clearly seen from this figure that the critical temperature reduces monotonically to zero by decreasing the single- 
ion anisotropy term D/\J\. For different values of the biaxial single-ion anisotropy, the critical temperature tend 
to zero for the ground-state value that is consistent with the ground-state phase boundary given by Eq. (27). 
While the uniaxial single-ion anisotropy term D forces spins to lie within x — y plane when D < Q, the E term 
tries to align them into the y — z plane. The biaxial single-ion anisotropy thus basically supports the magnetic 
ordering related to OP when D/\J\ < —1.5 and hence, it survives until more negative single-ion anisotropies 
D/\J\. On the other hand, the biaxial single-ion anisotropy additionally lowers the critical temperature of OP for 
the easy-axis single-ion anisotropy D > 0. It should be mentioned here that this behavior arises as a consequence 
of E term, which induces the quantum superposition between | + 1) and j — 1) spin states due to the nonzero 
quantum fluctuations arising from this term. Altogether, it could be concluded that the quantum reduction of 
the magnetization also lowers the critical temperature of OP owing to the quantum fluctuations closely associated 
with the E term. It should be also remarked that the critical temperatures of the mixed spin-1/2 and spin-1 
Ising model with only the uniaxial single-ion anisotropy term are recovered by neglecting the biaxial single-ion 
anisotropy, i.e. setting S = in Eq. (1). As a matter of fact, our results are in this particular case consistent 
with those reported on previously using other exact methods for the Bethe lattice models such as the non-linear 
mapping [24] , the cluster variational method in the pair approximation [5] or the method based on exact recursion 
relations [25, 26]. 

In order to provide a deeper insight into how the biaxial single-ion anisotropy affects the critical behavior, the 
critical temperature versus biaxial single-ion anisotropy dependence is shown in Fig. 4 for several values of the 
uniaxial single-ion anisotropy. Two different regions corresponding to OP and DP are separated by second- 
order transition lines. Evidently, it turns out that the biaxial single-ion anisotropy has major influence on the 
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critical behavior of the system. The critical temperature gradually decreases with increasing the biaxial single- 
ion anisotropy strength for D/\J\ > —1.0. It is quite clear that the suppression of critical temperature can be 
connected to the quantum fluctuations, which occur because of the effect of E term. In addition to this rather 
trivial finding, one also observes the interesting nonmonotonical dependences of the critical temperature. Namely, 
the critical temperature firstly increases and only then gradually decreases with a strength of the biaxial single-ion 
anisotropy for D/\J\ < —1.0. It should be pointed out that the spin-1 atoms are preferably thermally excited to 
the |0) state when D < 0, which means that they are preferably excited to the x — y plane. Since the biaxial 
single-ion anisotropy tries to align them into the y — z plane, it favors the spontaneous long-range order along 
2 axis in that it prefers the quantum superposition between | ± 1) spin states before the population of the |0) 
one. In agreement with the aforementioned arguments, the most interesting finding to emerge here is that there 
is a strong evidence for a spontaneous long-range order even under assumption of extraordinary strong (negative) 
single-ion anisotropies D/\J\ provided that the biaxial single-ion anisotropy is strong enough. The spontaneous 
long-range order arises in this rather peculiar case from the quantum fiuctuations caused by the biaxial single-ion 
anisotropy. When comparing the present exact results for critical temperatures with the ones for the similar 
mixed-spin Ising model on the honeycomb lattice [11], it can be concluded that the critical temperatures for the 
mixed-spin Ising model on the Bethe lattice are slightly greater than the ones for the honeycomb lattice. 
Now, let us study in detail the thermal dependences of magnetization that provide an independent check of the 
critical behavior. For easy reference, we will further use the extended Neel classification [7, 33] to refer to different 
types of temperature dependences of the total magnetization, which is schematically displayed in Fig. 5. The 
typical thermal variations of the sublattice and total single-site magnetizations are depicted in Fig. 6 for the 
uniaxial single-ion anisotropy D/\J\ — —1.0 and several values of the biaxial single-ion anisotropy E/\J\. As one 
can see from this figure, the sublattice magnetization uib varies very rapidly with E. If a strength of the biaxial 
single-ion anisotropy is rather large, the quantum reduction of the sublattice magnetization ms becomes stronger 
due to energetical favoring of the | — 1) spin state. The total magnetization then also changes from the standard 
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Figure 5. The extended Nccl classification for difTcrcnt types of temperature dependences of the total magnetization. Acronyms 
for particular cases are given by the circled letters. 
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Figure 6. (a) Temperature dependences of the sublatticc single-site magnetization for £>/| J| — —1.0 and several values of the 
biaxial single-ion anisotropy. (b) Temperature dependences of the total magnetization normalized per one site for 
D/\J\ — —1.0 and several values of the biaxial single-ion anisotropy. 



R-type dependence to the more interesting Q-type dependence as it can be seen in Fig. 6(b). 
As far as the biaxial single-ion anisotropy at the critical value E^/\J\ = \fTJ jl is considered, the sublattice 
magnetization mA fully compensates the sublattice magnetization ms in the ground-state as it can be seen from 
Fig. 2(b) and Fig. 7. In this particular case, the sublattice magnetization ttib becomes greater than the sublattice 
magnetization jtia, i.e. |ms| > ttia, for all nonzero temperatures below Tc. 

Finally, let us look more closely on another particular cases for different biaxial single-ion anisotropics close to 
the value E^. The total magnetization exhibits the P-type curves for E < a,s it is shown for E/\J\ — 2.57 and 
2.59. However, the sublattice magnetization \mA\ becomes greater than ms below Tc ii E > E^ is satisfied (see 
E/\J\ = 2.61 and 2.63). The more robust thermal excitations of the sublattice magnetization ms then result in the 
N-type dependence with one compensation point if the biaxial single-ion anisotropics are close enough but slightly 
higher than E^ (for instance E/\J\ — 2.61). It can be readily understood that the model exhibits one compensation 
point just in a very small restricted range of E/\J\ and therefore, the line of compensation temperatures is not 
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Figure 7. The various temperature dependences of the total magnetization normalized per one site when the strength of uniaxial 
single-ion anisotropy is fixed {D/\J\ — —1.0) and the biaxial single-ion anisotropy varies in the vicinity of 

shown for clarity in the finite-temperature phase diagrams. In the region, where the total magnetization exhibits 
one compensation point, the sublattice magnetization tua exceeds tub at lower temperatures, i.e. \mA\ > rns, 
while reverse is the case at higher temperatures, i.e. \mA\ < ran holds at higher temperatures (see E/\,J\ — 2.61). 
Last but not least, the R-type magnetization curves is obtained for the more stronger biaxial single-ion anisotropy 
such as E/\J\ = 2.63. 

4. Concluding remarks 

In this work, we have exactly studied the phase diagrams and magnetizations of the mixed spin-1/2 and spin-1 
Ising model on the Bethe lattice (q=3) by the use of two precise mapping transformations and exact recursion 
relations. The particular attention was focused on the effect of uniaxial and biaxial single-ion anisotropies acting 
on the spin-1 atoms. Exact results for the phase diagrams, total and sublattice magnetizations were obtained 
and discussed in detail. The obtained results show that the presence of a small amount of biaxial single-ion 
anisotropy has a very important influence on the magnetic properties. Beside this, we have found an exact 
evidence that the model exhibits within OP the remarkable quantum superposition of the | ± 1) spin states that 
leads to the quantum reduction of the magnetization whenever there is arbitrary but nonzero biaxial single-ion 
anisotropy. Macroscopically, this quantum superposition reduces the critical temperature of OP for the easy-axis 
uniaxial single-ion anisotropies {D > 0) as it also appreciably depresses the magnetization of spin-1 atoms from 
its saturation value even within the ground state (see Fig. 2(b)). The discussed reduction of critical temperature 
as well as magnetization is obviously of purely quantum origin, since it appears owing to the local quantum 
fluctuations arising from the biaxial single-ion anisotropy. On the other hand, the same quantum fluctuations 
can surprisingly cause an onset of the spontaneous long-range order related to OP for the extraordinary strong 
easy-plane anisotropies D < provided that there is strong enough biaxial single-ion anisotropy. Depending 
on a strength of the uniaxial and biaxial single-ion anisotropy parameters, the temperature variation of total 




Cesur Ekiz, Jozcf Strccka, Michal Jascur 



magnetization has been found to be either of Q-, R-, P-, L- or N-type. It should be also mentioned that the 
present results are in a good qualitative agreement with those obtained by using the same Hamiltonian on the 
honeycomb lattice [11]. 

Furthermore, it should be stressed that the present results are interesting both from the academic point of view 
(because of their exactness) as well as from the experimental viewpoint. Namely, the mixed-spin Ising models 
defined on the Bethe lattices can be useful in providing a deeper insight into a magnetism of dendrimeric organic 
molecules containing radicals as their magnetic sites [34-39]. In this respect, the lattice-statistical models with 
spin- 1/2 and spin-1 sites are the most useful ones as the most common dendrimeric organic magnets are based 
on radicals such as carbenes, polyethers, trityl and trityl-anionic radicals, perchlorotriphenylmethyl radical and 
so on, which are basic building blocks of either the spin- 1/2 or spin-1 sites. 

Finally, it is noteworthy that the exact solution of the spin- 1/2 Ising model on the triangular Husimi lattice has 
been obtained as another by-product of our calculations. Note that this interesting spin system represents an 
intermediate step of our exact procedure based on the star-triangle and triangle-star mapping transformations. 
According to Eqs. (15) and (16), exact results for the spin-1/2 Ising model on the triangular Husimi lattice with 
the nearest-neighbor interaction R can simply be found from the well-known exact results for the corresponding 
spin-1/2 Ising model on the Bethe lattice with the nearest-neighbor interaction K. For instance, the critical 
temperature of the spin-1/2 Ising model on the triangular Husimi lattice can be straightforwardly obtained by 
putting the critical temperature jScK = 2 In 3 of the spin-1/2 Ising model on the three-coordinated Bethe lattice 
into the mapping relation (15), which yields after elementary calculation 

/3.i? = lnQ^ ^M£^_L_ = 1.180225... (32) 

As could be expected, the critical temperature of the spin-1/2 Ising model on the triangular Husimi lattice with the 
effective coordination number g = 6 (i.e. the number of nearest-neighbors for each lattice site) lies in between the 
critical temperature of the spin-1/2 Ising model on the triangular lattice = = 0.910239 . . . and the critical 

temperature of the spin-1/2 Ising model on the six-coordinated Bethe lattice ^LaZk — = 1.233151 .... In 

this respect, the spin-1/2 Ising model on the triangular Husimi lattice slightly better approximates the spin-1/2 
Ising model on the triangular lattice. Finally, it is worthy to mention that exact procedure developed in this 
paper enables further remarkable extensions. For instance, the approach based on combining accurate mapping 
transformations with the exact recursion relations can be further generalized to account for the general spin-S* 
problem on the sublattice B, the next-nearest-neighbor interaction between the spin-1/2 sites from the sublattice 
A, the multispin exchange interaction between the spins from the sublattice A and B, etc. In this direction will 
continue our further work. 
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